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We present and use a technique for implementing in a fast way, and without any approximations, higher-order 
calculations of partonic cross sections into global analyses of parton distribution functions. The approach, which 
is set up in Mellin-moment space, is particularly suited for analyses of future data from polarized proton-proton 
collisions at RHIC, but not limited to this case. We also briefly discuss the effects of soft-gluon resummations on 
spin asymmetries. 



1. Introduction 

Current and future dedicated high-energy spin 
experiments are expected to vastly broaden our 
understanding of the nucleon spin structure. Pro- 
cesses in polarized lepton-nucleon and proton- 
proton scattering will be studied that will 
give access to, in particular, the spin-dependent 
gluon distribution Ag and the flavor-SU(2) break- 
ing in the polarized nucleon sea, Am — Ac?. 
Having available at some point in the near fu- 
ture spin data on various different reactions, one 
needs to tackle the question of how to determine 
such polarized parton densities from the measure- 
ments. This is not at all a new problem: in 
the unpolarized case, several groups perform such 
"global analyses" of the plethora of data available 
there |||]]. The strategy is in principle clear: an 
ansatz for the parton distributions at some initial 
scale fio, given in terms of appropriate functional 
forms with a set of free parameters, is evolved to 
a scale /if relevant for a certain data point for a 
certain cross section. Then the parton densities 
at scale [If are used to compute the theoretical 
prediction for the cross section, and a x 2 value is 
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assigned that represents the quality of the com- 
parison to the experimental point. This is done 
for all data points to be included in the analysis, 
and subsequently the parameters in the ansatz for 
the parton distribution functions are varied, until 
eventually a minimum in \ 2 is reached. 

In practice, this approach is difficult if the par- 
tonic scattering is treated beyond the lowest or- 
der of perturbation theory (PT) . The numerical 
evaluation of the hadronic cross section at higher 
orders is usually a rather time-consuming proce- 
dure as it often requires several tedious numerical 
integrations, not only for the convolutions with 
the parton densities, but also for the phase space 
integrations in the partonic cross section. The fit- 
ting procedure outlined above, on the other hand, 
usually requires thousands of computations of the 
cross section for any given data point, and so the 
computing time required for a fit easily becomes 
excessive even on modern workstations. 

In the unpolarized case, a way to get around 
this problem is based on the fact that the par- 
ton densities are already known here rather accu- 
rately HQ] . As a consequence, the theory answer 
for a certain cross section is expected to change 
in a very predictable way when going from the 
lowest-order Born level to the first-order approx- 
imation. It is then possible to pre-calculate a set 
of correction factors Ki (i running over the data 
points) , and to simply multiply them in each step 
of the fitting procedure to the lowest-order ap- 
proximation for the cross section, the latter being 
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usually much faster to evaluate than that involv- 
ing higher order terms. In the polarized case, 
however, it is in general not at all clear whether 
such a strategy will work. Here, the parton densi- 
ties are known with much less accuracy so far. It 
is therefore not possible to pre-calculate higher- 
order correction factors that one would be able to 
keep fixed throughout the fit, while using "fast" 
lowest-order expressions for the partonic cross 
sections. For instance, the spin-dependent par- 
ton distributions, as well as the polarized partonic 
cross sections, may have zeros in the kinematical 
regions of interest, near which the predictions at 
lowest order and the next order will show marked 
differences. 

Our goal is to find a way of implementing 
efficiently, and without approximations, higher- 
order expressions for any hadronic cross sec- 
tion into the fitting procedure. As was recently 
shown ||, this can be achieved in a very sim- 
ple and straightforward way by going to Mellin-n 
moment space. A related technique was first dis- 
cussed in HQ. 

In the following we will outline our method 
and apply it to an example, the production of 
prompt photons at high transverse momentum in 
polarized pp collisions at RHIC. The sensitivity 
of this reaction to the gluon distribution via the 
LO Compton subprocess is the main reason why 
this process will be the flagship measurement of 
Ag at RHIC [Q. As a first case study for future 
global analyses we also carry out a "toy" analy- 
sis of polarized deep-inelastic scattering (DIS) Q 
and prompt photon data projected for RHIC. 

At the end of this paper, we digress from our 
main topic and briefly discuss the effects of soft- 
gluon resummation on spin asymmetries. 

2. Hadronic Cross Sections and the Mellin 
Moment Technique 

A general spin-dependent cross section in longi- 
tudinally polarized pp collisions, differential in a 
certain observable O and integrated over experi- 
mental bins in other kinematical variables T, can 
be written as 
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where the arguments (-1 — h) and (H — ) in the first 
line refer to the helicities of the incoming protons. 
The Afi are the spin-dependent parton distribu- 
tions, defined as 



Afi(x,fi F ) = ff(x,n F ) -f t (x,fj, F ) 



(2) 



where /. (fl~) denotes the number density of a 
parton-type fi with hclicity '+' (' — ') in a pro- 
ton with positive hclicity, carrying the fraction x 
of the proton's momentum. Note that for some 
observables, a fragmentation function may be ad- 
ditionally present in Eq. (|l|). 

The scale fi F is the factorization scale for initial 
state collinear singularities and reflects the cer- 
tain amount of arbitrariness in the separation of 
short-distance and long-distance physics embod- 
ied in Eq. ([!]). The other scale, [jlr, in Eq. ([!]) 
is the renormalization scale, introduced in the 
procedure of renormalizing the strong coupling 
constant. Finally, the sum in Eq. ([!]) is over all 
contributing partonic channels a + b — > H + X, 
with dAa^ b the associated spin-dependent par- 
tonic cross section, defined in complete analogy 
with the first line of Eq. ([!]) . They are perturba- 
tive and hence are expanded in powers of a s . 

The Mellin moments of the polarized parton 
distribution functions are defined as 

A/f( M )= / dxx^AMx^) ■ (3) 
Jo 

The parton distributions in Bjorken-ir. space are 
recovered from the moments by an inverse Mellin 
transform: 



Af i {x,n) = ±-. [ dnx- n Af^) 



(4) 



where C n denotes a contour in the complex n 
plane that has an imaginary part ranging from 
— oo to oo and that intersects the real axis to the 
right of the rightmost poles of the A/™(^i). 

The crucial, but simple, step in applying mo- 
ment techniques to Eq. (fy) is to express the 
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Afi(xi, hf) by their Mellin inverses in Eq. (|J) 
(7||| . One subsequently interchanges integrations 
and arrives at 
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One can now pre-calculate the quantities 
Aa^ b (n,m,0, (j,r, (if), which do not depend at 
all on the parton distribution functions, prior to 
the fit for a specific set of the two Mellin vari- 
ables n and m, for each contributing subprocess 
and in each experimental bin. Effectively, one 
has to compute the cross sections with complex 
"dummy" parton distribution functions x~ n x b m . 
All the tedious and time-consuming integrations 
are already dealt with in the calculation of the 

Aa^ b (n,m,0,fj, R ,n F )- 

The double inverse Mellin transformation 
which finally links the parton distributions 
with the pre-calculated Aa^ b (n,m,0, fin, fip) of 
course still needs to be performed in each step 
of the fitting procedure. However, the integra- 
tions over n and m in Eq. (^) are extremely 
fast to perform by choosing the values for n, m 
in Aa^ b (n, m, O, fiR, ^lf) on the contours C n , C m 
simply as the supports for a Gaussian integra- 
tion. The point here is that the integrand in n 
and m falls off very rapidly as \n\ and \m\ in- 
crease along the contour, for two reasons: first, 
each parton distribution function is expected to 



fall off at least as a power (1 



at large x, 



which in moment space converts into a fall-off of 
~ l/n 4 or higher. Secondly, we may choose con- 
tours in moment space that are bent by an angle 
a — it/2 with respect to the vertical direction; a 
possible choice is shown in Fig. [j]. Then, for large 
\n\ and \m\, n and m will acquire large negative 
real parts, so that (x a )~ n and (x b )~ m decrease 
exponentially along the respective contours. This 



helps for the numerical convergence of the calcu- 
lation of the Aa^ b (n, m,0, fin, /if) and also gives 
them a rapid fall-off at large arguments. 




Figure f. Contours in complex Mellin-n, m spaces 
for the calculation of the double Mellin inverse in 
Eq. (||). r n and r m denote the rightmost poles of 
the integrand in n and m, respectively, and the Cj 
the intersections with the real axis. 



3. Example: Prompt Photons at RHIC 

In Rcf. H we have presented two practical 
applications of the Mellin technique: the semi- 
inclusive production of hadrons in polarized DIS, 
and the production of prompt photons in polar- 
ized pp collisions at RHIC. Here we will briefly 
recall our findings for the latter process. To 
be specific, we consider the transverse momen- 
tum (j>t) distribution of the prompt photon, in- 
tegrated over a certain experimental bin in its 
pseudorapidity r\. Thus, we have "O = pt" and 
"T = 77" in Eq. (§). 

For prompt photon production, the lowest- 
order partonic reactions are q + q — * 7 + 5 and 
q + g —>■ 7 + q, the latter channel being sensitive 
to the polarized gluon distribution. The next- 
to-leading order (NLO) corrections are also avail- 
able U and will be used in our case study. 

We use VS = 200 GeV and consider five values 
of pt which will be experimentally accessible at 
RHIC, p T = [12.5,17.5,22.5,27.5,32.5] GeV. We 
average over \rj\ < 0.35 in pseudorapidity and im- 
pose an isolation cut on the photon 0] , for which 
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we choose the criterion proposed in |lCf | with pa- 
rameters R — 0.4, e = 1. A positive feature of this 
isolation criterion is the absence of a fragmenta- 
tion contribution to prompt photon production. 
We choose the renormalization and factorization 
scales iir — (j,f = Pt- 

Our first goal here is to show that the method 
based on Eq. (j^) actually works and correctly re- 
produces the result obtained within the direct, 
but "slow", calculation via Eq. d). Also, we 
need to establish an optimal size of the grids 
that yields excellent accuracy but is still calcu- 
lable in, say, a few hours of CPU time on a stan- 
dard workstation. Fig. compares the results 
based on Eq. (||), referred to as "Mellin tech- 
nique" , to those of Eq. (Q) , for various sizes of 
the grid in n,m. Here we have used the polar- 
ized parton densities of |ll[] ("standard" set). 
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Figure 2. Comparison of the results based on the 
Mellin-technique in Eq. ([l]) to those of Eq. (||) for 
various sizes of the grid in n,m. 

For a more detailed comparison, we split up the 
contributions to the NLO prompt photon cross 
section into three parts, associated with the reac- 
tions q + q — > 7 + X and q + g — > 7 + X that are 
already present at Born level, and all other pro- 
cesses that arise only at NLO. One notices that 
in each case already a grid size of 64 x 64 values 
yields excellent accuracy. Even a 56 x 56 grid is 
acceptable apart from a minor deviation occur- 
ring for qq scattering in the vicinity of a zero in 
the partonic cross section. 

The crucial asset of the Mellin method is 
the speed at which one can calculate the 



full hadronic cross section, once the grids 
Aa2i,(n,m,pT, [ir, hf) have been pre-calculated. 
For the 64 x 64 grid, we found that 1000 evalu- 
ations of the full NLO prompt photon cross sec- 
tion take only about 10-15 seconds on a standard 
workstation. Note that this number includes the 
evolution (in moment space) of the parton dis- 
tributions from their input scale to the scale pr 
relevant to this case. Clearly, an implementation 
into a full parton density fitting procedure is now 
readily possible. 

To give an example, we finally perform a "toy" 
global analysis of the available data on polar- 
ized deep-inelastic scattering || and of fictitious 
data on prompt photon production at RHIC 
which we project by simply calculating A~' LL to 
NLO using the sets of polarized and unpolarized 
parton distributions of |Tl]] and respectively. 
For an estimate of the anticipated la errors on 
the "data" for A 1 LL , we use the numbers reported 
in ||. We subsequently apply a random Gaus- 
sian shift of the pseudo-data, allowing them to 
vary within la. The "data" , as well as the un- 
derlying theoretical calculation of A 7 LL based on 
the spin-dependent parton densities of |ll|] (solid 
line), are shown in Fig. [|(a). 

Next, we perform a large number of fits to 
the full, DIS plus projected prompt photon, "data 
set" . We simultaneously fit all polarized parton 
densities, choosing the distributions of [|ll| as the 
input for the Aq, Aq, but using randomly chosen 
values for the parameters in the ansatz for the po- 
larized gluon distribution at the input scale fj,Q. 
Regarding the details of the evolution, we stay 
within the setup of |Q, but we choose a more 
flexible ansatz for Ag: 

Ag{x, ti )=Nx a (l- x) 13 (1 + 7 x) g(x, fi ) , (6) 

which also allows for a zero. g(x, /iq) is the unpo- 
larized gluon density Q at the input scale of (ll]] . 
Ideally, thanks to the strong sensitivity of the 
prompt photon reaction to Ag, the gluon density 
in each fit should return close to the function we 
assumed when calculating the fictitious prompt 
photon "data" , in the region of x probed by the 
data. Indeed, as shown in Fig. ||(b), this hap- 
pens. The shaded band illustrates the deviations 
of the gluon densities obtained from the global 
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fits to the "reference Ag" [jll| used in generat- 
ing the pseudo-data. It should be stressed that 
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As is expected, the gluon density is rather 
tightly constrained in the ^-region dominantly 
probed by the prompt photon data. This is true 
in particular atiw 0.15, as a result of the most 
precise "data point" for A\ L at pt = 12.5 GeV. 
Fig. ||(b) shows also two extreme gluon densi- 
ties with first moments Ag 1 (fi ) = ±0.8 (dot- 
ted lines), which are both in perfect agreement 
with all presently available DIS data. It should be 
noted that future measurements of A 1 LL at RHIC 
at \/~S — 500 GeV and for similar pt values would 
further constrain Ag in the x-region between 0.05 
and 0.1. Although our analysis still contains a 
certain bias by choosing only the framework of 
]Tl| ] for the fits as well as by our choice of what 
X 2 values are still tolerable, it clearly outlines the 
potential and importance of upcoming measure- 
ments of A\ L at RHIC for determining Ag. 
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Figure 3. (a): Generated pseudo data for A 1 LL 
based on a calculation using the spin-dependent par- 
ton densities of [|ll) (solid line). The shaded band 
corresponds to the results of a large number of com- 
bined fits to DIS and AJ^ L data (see text), (b): 
Variations of the polarized gluon densities obtained 
in the combined fits, with respect to Ag of [jljj . 

only those fits are admitted to the band that give 
a good simultaneous description of the DIS and 
A~l L data. Here we have tolerated a maximum 
increase of the total \ 2 by up to four units from 
its minimum value. The shaded area in Fig. ||(a) 
shows the variations in Al T . 



4. Spin asymmetries and resummation 

In the unpolarized pattern of disagree- 

ment between theoretical predictions and exper- 
imental data for prompt photon production has 
been observed in recent years not globally 
curable by 'fine-tuning' the gluon density Jl3| . 
The main problems reside in the fixed-target re- 
gion, where NLO theory dramatically underpre- 
dicts some data sets. At collider energies, as rel- 
evant to RHIC, there is less reason for concern. 

In view of this, various improvements of the 
theoretical framework have been developed. One 
of them resorts to applying 'threshold' resum- 
mation to the prompt photon cross section Jl4[ , 
which organizes to all orders in a s large loga- 
rithmic corrections to partonic hard scattering 
associated with emission of soft gluons. As the 
partonic c.m.s. energy s approaches its mini- 
mum value at s = ip^, corresponding to 'par- 
tonic threshold' when the initial partons have just 
enough energy to produce the photon and the re- 
coiling jet, the phase space available for gluon 
bremsstrahlung vanishes, resulting in corrections 
to the partonic cross section da/dpr as large as 
c^m 2fe (l - 4p^/s)a Bmn at fc-th order in PT. 
Threshold resummation jL5|,[l4j organizes this sin- 
gular behavior of dcr/dpT to all orders. It is again 
carried out in Mellin-n moment space, where the 
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Figure 4. Spin asymmetry All at LO, NLO, and 
including NLL threshold resummation. 

logarithms have the form a* ln 2fe (n) a Born (n) . Its 
application is particularly interesting in the fixed- 
target regime, since here the highest xt are at- 
tained in the data and the discrepancy between 
data and theory is largest. 

In phenomenological applications of 
threshold resummation, one finds a significant, 
albeit not sufficient, enhancement of the theory 
prediction in the fixed-target regime at large val- 
ues of pt I VS, accompanied by a dramatic re- 
duction of scale dependence [^6[0. Thanks to 
the universal structure of soft-gluon emission, it 
is straightforward to apply threshold resumma- 
tion to the polarized cross section. Fig. |] shows 
the resulting effects on the spin asymmetry All-, 
for a 'toy' example that assumes a fictitious po- 
larized set-up of the E706 experiment. Details 
are as in |l6| ]. Even though Fig. ^ does not di- 
rectly refer to the case of RHIC, it is good news 
that resummation effects cancel to a large extent 
in All for our present example. 

A similar cancellation of higher-order effects as- 
sociated with soft-gluon emission is found for the 
parity-violating longitudinal single-spin asymme- 
try for jet production in pp collisions, which has 
been proposed Jl8|,|| as a candidate indicator for 
physics beyond the Standard Model and is also 
measurable at RHIC. 
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